# autoplot(prcomp(t(mcp_zscore_matrix)), data = person_metadata_df, colour = 'person', main = "MCPCounter")
# autoplot(prcomp(t(mcp_zscore_matrix)), data = person_metadata_df, colour = 'cell_type', main = "MCPCounter")
# autoplot(prcomp(t(cs_zscore_matrix)), data = person_metadata_df, colour = 'person', main = "Cibersort")
# autoplot(prcomp(t(cs_zscore_matrix)), data = person_metadata_df, colour = 'cell_type', main = "Cibersort")
pheatmap(
mcp_zscore_matrix,
main = "MCPCounter genes",
annotation_row = mcp_heatmap_row_df,
annotation_col = heatmap_col_df,
scale = "none")

pheatmap(
cs_zscore_matrix,
main = "Cibersort genes",
annotation_col = heatmap_col_df,
scale = "none")

# pheatmap(
# mcp_mean_zscore_matrix,
# main = "MCPCounter genes, collapsed by sample/celltype",
# annotation_row = mcp_mean_heatmap_row_df,
# annotation_col = mean_heatmap_col_df,
# scale = "none")
# pheatmap(
# cs_mean_zscore_matrix,
# main = "Cibersort genes, collapsed by sample/celltype",
# annotation_col = mean_heatmap_col_df,
# scale = "none")
# kable(mcp_result_df)
# kable(cs_result_df)